Terahertz radiation due to random grating coupled surface plasmon polaritons 
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We report on terahertz (THz) radiation under electrical pumping from a degenerate semiconduc- 
tor possessing an electron accumulation layer. In InN, the random grating formed by topographical 
defects provides high-efficiency coupling of surface plasmon polaritons supported by the accumula- 
tion layer to the THz emission. The principal emission band occupies the 2-6 THz spectral range. 
We establish a link between the shape of emission spectra and the structural factor of the random 
grating and show that the change of slope of power dependencies is characteristic for temperature- 
dependent plasmonic mechanisms. The super-linear rise of a THz emission intensity on applied 
electric power provides advantage of such materials in emission yield. 

PACS numbers: 73.20.Mf 



The surface plasma waves confined to the interface be- 
tween a conductor (metal or degenerate semiconductor) 
and a dielectric attract a lot of attention for terahertz 
(THz) frequency range application [l|, Q. The character- 
istic frequency of a surface plasmon at the sharp interface 
of the semi-infinite conductor with vacuum reads as 



= 47r7VeVm*(£: + l), 



(1) 



where TV, e, and m* are the bulk concentration, charge, 
and effective mass of carriers, respectively; e is the dielec- 
tric permittivity of the conductor. This frequency is in 
the visible range for the metals, being in the far infrared 
and even in the terahertz ranges for the degenerate semi- 
conductors The distinct difference between plas- 
mon performances in these two types of materials is that 
the strong interband absorption in the visible range ham- 
pers collective electron excitations in the metals, while 
in the semiconductors the plasmon frequency is below 
an absorption edge. Thus, long-range surface plasmon 
polariton modes can be supported. 

For strictly two-dimensional plasmons, the frequency 
depends on the plasmon wave vector k as cj(/c) oc 
(sf. (This relation is more complicated for plasmons in 
the films of a finite thickness (Jlsl.) It determines why 
the surface plasmons cannot radiate to ambient directly. 
Their wave vector is always higher than that of light 
ko = cj/c at the same frequency. Coupling of the plasmon 
with an electromagnetic wave occurs by means of the mo- 
mentum exchange either at the disturbances induced by 
plasma wave instability JlJ or at the periodical modula- 
tions of surface profile ^]^. For k ^ k^., the radiating 
states have the momentum k ^ 27rm/a (m = 1,2,...), 
where a is a grating period. In the case of the ideal 
grating, the k intervals where the emission yield takes 
place are inevitably narrow and the integral emission in- 
tensity is not sufficient sometimes for reliable registration 
by available detectors. To circumvent this obstacle, more 
complicated (multiple-sections) gratings have been pro- 
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FIG. 1: Dispersion dependencies (semi- logarithmic scale) cal- 
culated for the different thicknesses d of plasmon supporting 
layers shown together with the light wave dispersion inside 
InN. The crossing of vertical dashed lines with plasmon dis- 
persion curves indicates the d thickness range suitable for the 
first order coupling {k ^ 27r/a) for a — 0.1-^1 iim. 



posed [ilTll. 

In this letter, we report on THz radiation from a de- 
generate semiconductor possessing a carrier accumula- 
tion layer which can support surface plasmon polaritons. 
We demonstrate advantages of the random grating to 
provide effective coupling between these plasmons and 
the THz emission in a wide spectral region extended for 
several THz. The genesis of the THz emission is con- 
sidered to establish the interplay of the distribution of 
grating wave vectors and a spectra shape. It is eluci- 
dated how the thermal population of plasmonic states 
controls the spectra and power dependencies. 

The study has been done using InN epilayers which 
have an intrinsic electron concentration of 10^^ — 10^^ 
cm~^, that guarantees uOs within a terahertz range. InN 
has a tendency to the formation of various structural 
defects, such as pores, pits, and trenches, frequently 
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decorated by metallic nanoparticles. The average dis- 
tance between these defects is about the grain size ('^1 
/im); therefore, a certain periodicity exists in their to- 
pographic position. The 1-3 /im thick InN layers were 
grown by plasma-assisted molecular beam epitaxy on c- 
sapphire atop of GaN buffers at substrate temperature 
Ts = (450 - 550) °C. For the THz electroluminescence 
measurements, contacts are formed on the top surface 
of samples. The THz spectra were measured at 8 K us- 
ing a step-scan Fourier spectrometer with excitation by 
the series of packets of 15 V rectangular pulses of 10 jas 
duration (71.5 Hz repetition rate). The dependencies of 
integral THz emission intensity on the electric field were 
measured at 4.2 K using 2-/is-long pulses by a Ge:Ga 
photodetector with 2.4-5 THz spectral sensitivity ran^e. 
Other experimental details are described elsewhere |12| . 

Calculations based on numerical solving of the Maxwell 
equations for multiple layers structures have been used 
to model the full set of plasmonic and waveguide modes 
supported in the studied samples. They have shown that 
the range of plasmonic layer thicknesses corresponding 
to the experimentally observed emission frequencies of 
several THz is rather narrow. For the sake of demonstra- 
tion, we present in Fig. 1 the plasmon dispersion curves 
calculated using the three-layers model as reported pre- 
viously [1, [HI. In these calculations, the metal-like 
layer described by the Drude model {N = 10^^ cm~^, 
m* = 0.13mo, damping parameter of 1 meV) is placed 
between the vacuum and the rest part of the InN epilayer, 
which is considered as a dielectric. The static permittiv- 
ity of InN was assumed equal to 14 everywhere. The 
results of the calculations explicitly display that to pro- 
vide the coupling by the typical InN defects, the layer 
thickness d should be small enough, < 0.1 /im. This 
estimation is consistent with similar limitation on the 
d value needed to support the long propagating surface 
modes 0, El • A good candidate for such a thin plasmon- 
supporting layer is an electronjLCCumulation layer on a 
top of a semiconductor film |13j. The presence of the 
layer in InN has been confirmed by the spectroscopy of 
electron energy losses induced by the surface plasmons 

Here, we demonstrate other observation related pre- 
sumably to the surface plasmon-polaritons, which is ob- 
tained by scanning electron microscopy (SEM). Before 
to its description, it is worth reminding that the elec- 
tron beam excitation is a very efficient way to create 
the surface plasmons [lH, and that this excitation is 
followed by the generation of secondary electrons (SEs) 
due to non-radiative decay of plasmons into single elec- 

The SEM studies were done 



tron excitations 
at low temperature to minimize the plasmon damping. 
The electron beam impinges normally to a cleaved facet 
of a layer. This configuration corresponds to the so- 
called aloof scanning mode, when the electrons moving 
along nontouching trajectory above the sample surface 
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FIG. 2: Facet SEM image of an InN epilayer recorded at 5 
K using a SE2 detector of a SEM analytical microscope (20 
keV electron beam) . The inset presents the image taken with 
a sample detuned slightly to show discretion in a cloud base. 



can stimulate the surface plasmons in conductors 16|, [l^ • 
The SEM images of the InN epilayers exhibit intricate 
features - clouds above the sample surface, which have a 
specific semispherical shape (Fig. 2). Their sizes increase 
with exposure to the electron beam. With slight declin- 
ing a sample, it was possible to resolve discrete structure 
in their base, whose step is approximately equal to the 
average distance between the topographic defects. The 
plasmon excitation in the aloof mode implies that the 
defects scatter the electrons with consecutive charge ac- 
cumulation. We suppose that the secondary electrons 
display an electric field around the charged defects in the 
same way as the iron filings depict the magnetic field 
lines. Since SEs can escape from rather small depth of 
several nanometers only, the registration of such clouds 
confirms implicitly the propagation of the plasma waves 
in a very thin electron accumulation layer on a top of a 
semiconductor. 

To consider the genesis the THz emission from the 
InN layer under electrical pumping, we assume first that 
each harmonic of the random grating translates the cor- 
responding plasmons into a light cone at the frequency 
governed by both the plasmon dispersion and the par- 
ticular wave vector. To find the structure factor of the 
random grating, S{k), the Fourier transformation of a 
SEM image was performed [Fig. 3 (a, b)]. Then, it is 
quite reasonable to suggest that the intensity of the radi- 
ation at certain energy is proportional to the amplitude 
of respective harmonic. The influence of the plasmon 
dispersion, calculated as described above, is taken into 
account by multiplying on the plasmon density of states 
in this point. Such an approach permits us to simulate 
directly an emission spectrum, neglecting higher order 
scattering. Figure 3 (c) demonstrates the good agree- 
ment between the computed and experimental spectra, 
in spite of a rather complicated shape provided by the 
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FIG. 3: (color online) Simulation of THz emission spectra. 

(a) SEM image of a surface of an InN epilayer (sample A). 

(b) Image of the structure factor of random grating obtained 
by the Fourier transformation of the SEM image. Color scale 
gives the harmonic amplitudes. Dark spot in the center is due 
to removing the unnecessary harmonics of lowest orders, (c) 
Comparison of experimental and calculated [d —7 nm, other 
parameters are given in the text) spectra. 



random grating. 

In accord with the energy and momentum conserva- 
tion requirements, the higher energy electrons populate 
plasmonic states with lower k- vectors [18|. In our calcu- 
lations, we ignore such a nonequilibrium plasmon popu- 
lation that provides, likely, some underestimation of the 
intensity of the lower-energy peaks in the calculated spec- 
trum as compared with the experimental one. The other 
difference in the spectra concerns the higher energy cut- 
off of the main emission band, which is at higher energy 



in the simulated spectrum. The difference arises out dis- 
regarding the finite population of the plasmonic states. 
At electrical pumping, the filling level is controlled by 
the electron temperature T, which can markedly exceed 
the lattice temperature [l, [l9| . Its value is roughly pro- 
portional to the electrical power P^i. Our spectral mea- 
surements were done using a spared regime, i.e. with the 
small Pez; thus, the actual filling level was lower than the 
energy corresponding to the largest /c^^^ in the grating 
wave vector distribution. 

The filling effect influences also the dependencies of 
emission intensity on electric power, namely, their char- 
acteristic slopes. To describe that, we suppose the plas- 
mons to be thermalized with temperature T, e.g., via 
the electron scattering by phonons, impurities, and other 
defects [IqI, |2Q|. The emitting power W{T) can be esti- 
mated by integration over all plasmon states with the 
energy E{k)^ occupation of these states, and the struc- 
tural factor of the grating S{k)'. 



W(T) = J kdkS{k)E{k)/{eyip{E{k)/T) 



1). (2) 



Here, the temperature T is given in the energy units. 
Two following limiting cases are considered: i) homoge- 
neous emission from all possible plasmon states assuming 
5'(/c)=const; ii) emission from the grating-coupled states 
with S{k) derived from the SEM image (see Fig. 3). 

The performed calculations show that the mechanism 
based on the uniform scattering dominates when T < 
£^(/c^*^), i.e. when the states corresponding to the min- 
imal grating wave vector are not populated yet [Fig. 4 
(a)]. Consequently, the slope in the resulting emission 
power dependency must have the exponent /3 ~ 5. (This 
law appears due to the superposition of the Planck 
distribution function and the plasmon dispersion.) With 
the temperature rise, the mechanism of the emission 
changes - now the most of the emission comes from the 
plasmons interacting with the random grating. When 
the temperature approaches the higher energy cut-off 
T - £;(A;^^^), the /3 - 1.3 is expected. The predicted 
slopes and their change do be observed in the experimen- 
tal dependencies of the radiated intensity on the electric 
power [Fig. 4 (b)]. Note that the distinctiveness of the 
kink and the corresponding P^i value are specific for each 
sample, being dependent on the grating pattern and plas- 
mon dispersion. 

For the sake of clarity, the presented consideration is 
markedly simplified. Nevertheless, it displays the basic 
features of the plasmonic emission from nanostructured 
degenerate semiconductors. In particular, the super- 
linear rise of the THz intensity from InN with increasing 
electric power (/3 > 1.3) prevails over sub-linear laws in- 
herent to other processes in semiconductor structures at 
high electric field [l2|, [2lj. As a result, the THz emission 
with the density as high as 30 /iW/cm^ can be collected 
from InN epilayers. The 2-6 THz spectral range occupied 
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FIG. 4: (color online) (a) Dependencies of the emission in- 
tensity on temperature (given in energy units) calculated for 
the homogeneous filling of plasmon states and for the filling 
of grating-coupled states. The vertical lines mark the interval 
corresponding to the occupation of the grating-coupled states 
in the sample A. At low temperature, when the emission pro- 
vided by the homogeneous filling dominates, the exponent /3 
is about 5; at the higher-energy boundary, it drops to ^1.3. 
(b) Experimental dependencies of the emission intensity on 
the electrical power of representative samples exhibiting the 
similar slopes and the kink between them. The lines /3 = 5 
and /3 = 1.3 are given to guide the eye. 

by the emission band overlaps the range of the spectral 
sensitivity of a standard detector (2.4-5 THz). Thus, full 
ability of such detectors can be exploited using an emitter 
with random grating. 

In conclusion, our studies of InN have exhibited that 
the degenerate semiconductors possessing an electron ac- 
cumulation layer can produce terahertz radiation related 
to surface plasmon polaritons. An electron temperature 
increase, provided by applying electric field, populates 
various electronic states. In the plasmon momentum dis- 
tribution, equilibrium is broken down by the presence 
of the random grating formed by topographic defects. 
Thus, THz emission of appreciable intensity appears in 
the spectral range determined by both plasmon disper- 
sion and statistics of grating wave vectors. We believe 
that similar nanostructured materials can be promising 
for various THz application. 
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